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ABSTRACT 

We present a measurement of the angular power spectrum of the cosmic microwave background 
(CMB) using data from the South Pole Telescope (SPT). The data consist of 790 square degrees 
of sky observed at 150 GHz during 2008 and 2009. Here we present the power spectrum over the 
multipole range 650 < i < 3000, where it is dominated by primary CMB anisotropy. We combine 
this power spectrum with the power spectra from the seven-year Wilkinson Microwave Anisotropy 
Probe (WMAP) data release to constrain cosmological models. We find that the SPT and WMAP 
data are consistent with each other and, when combined, are well fit by a spatially flat, ACDM 
cosmological model. The SPT+WMAP constraint on the spectral index of scalar fluctuations is 
Us = 0.9663 ± 0.0112. We detect, at ~5cr significance, the effect of gravitational lensing on the CMB 
power spectrum, and find its amplitude to be consistent with the ACDM cosmological model. We 
explore a number of extensions beyond the ACDM model. Each extension is tested independently, 
although there are degeneracies between some of the extension parameters. We constrain the tensor- 
to-scalar ratio to be r < 0.21 (95% CL) and constrain the running of the scalar spectral index to be 
dns/d\nk = —0.024 ± 0.013. We strongly detect the effects of primordial helium and neutrinos on 
the CMB; a model without helium is rejected at 7.7a, while a model without neutrinos is rejected 
at 7.5(7. The primordial helium abundance is measured to be 1^ = 0.296 ± 0.030, and the effective 
number of relativistic species is measured to be TVcff — 3.85 ± 0.62. The constraints on these models 
are strengthened when the CMB data are combined with measurements of the Hubble constant and 
the baryon acoustic oscillation feature. Notable improvements include Ug = 0.9668 ± 0.0093, r < 0.17 
(95% CL), and N^s = 3.86±0.42. The SPT-I-WMAP data show a mild preference for low power in the 
CMB damping tail, and while this preference may be accommodated by models that have a negative 
spectral running, a high primordial helium abundance, or a high effective number of relativistic species, 
such models are disfavored by the abundance of low-redshift galaxy clusters. 

Subject headings: cosmology - cosmologyxosmic microwave background - cosmology: observations - 
large-scale structure of universe 
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1. INTRODUCTION 

Measurements of anisotropy in the temperature of 
the cosmic microwave background (CMB) are among 
the most informative and robust probes of cosmology. 
The acoustic oscillations of the primordial plasma have 
been measured on degree scales {£ < 500) with cosmic- 
variance-limited precision by the Wilkinson Mic rowave 



Anisotropy Probe (WMAP) ( [Larson et a L| [ "2011 [), yield- 
ing a wealth of cosmological intormation ( [Komatsu et al. 



2011 1. On much smaller scales, £ > 3000, the millimeter 



wave anisotropy is dominated by secondary anisotropics 
from the Sunyaev-Zel'dovich (SZ) effects and by emis- 
sion from foreground galaxies. The thermal SZ effect 
arises from the scattering of CMB photons off the hot 



gas in gra.vitatio nally collapsed structures (Sunyaev & 



Zel'dovich 197^, and thereby encodes information on 



the amplitude of matter fluctua t ions a t intermediate red- 
shifts. Recently, Lueker et al. (2010) reported the first 



statistical measurement of the SZ effect using multi- 
frequency South Pole Telescope (SPT) data. This was 
followed by a measurement u sing data from th e Atacama 
Cosmology Telescope (ACT, Das et al. 2011a; Dmikley 
et al. 2010l) and a n improved SP'l' measurement ( ,Shi-| 
rokott et al. 2011). The angular power spectrum of 
millimeter- wave emission from high-redshift, dusty, star- 
forming galaxies has also bee n characteri ze d by SPT, 
ACT, and Planck (jHall et al. 2010; Dunklcy et al.||2010l 



Shirokoff et al. 2011, Planck Collaboration 20 
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(Jn intermediate scales, 500 < t < 3000, the primary 
CMB anisotropy is the dominant source of millimeter- 
wave anisotropy, but its power is falling exponentially 
with decreasing angular scale. The reduction in CMB 
power is due to the diffusion of photons in the primor- 
dial plasni a and is often referred to as Silk damping 
(iSilkl 19681). This "damping tail" of the primary CMB 



anisotropy has been measured by a number of experi 
ments, notably the Arc minute Cosmology Bo lometer Ar 
ray Receiver (ACBAR,lReic hardt et al.|2009|), QUEST a t 
DASI (QUaD, Brown et al. 2009| |}''riedman et al.||2009[ ), 
and ACT (Das et al. 2011a i ) . 

Measurements ot tne CMB damping tail, in conjunc- 
tion with WMAP's measurements of the degree-scale 
CMB anisotropy, provide a powerful probe of early- 
universe physics. The damping tail measurements sig- 
nificantly increase the angular dynamic range of CMB 
measurements and thereby improve the constraints on 
inflationary parameters such as the scalar spectral in- 
dex and the amplitude of tensor fluctuations. Measure- 
ments of the angular scale of the damping can constrain 
the primordial helium abundance and the effective num- 
ber of relativistic particle species during the radiation- 
dominated era. Finally, the damping tail is altered at the 
few-percent level by gravitational lensing of the CMB, 
and is therefore sensitive to the matter fluctuations at 
intermediate redshifts. 

The work presented here is a measurement of the CMB 
damping tail using data from the SPT. The data were 
taken at 150 GHz during 2008 and 2009 and cover ap- 
proximately 790 square degrees of sky. This is approx- 
imately four times the area used in the preceding SPT 
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power spectrum result, Shirokoff et al. (2011). The new 
power spectrum spans the multipole range 650 < £ < 
3000 (angular scales of approximately A' < 9 < 16') and 
is dominated by primary CMB temperature anisotropy. 

The paper is organized as follows. We describe the 
SPT, the observations used in this analysis, and the 
pipeline used to process the raw data into calibrated 
maps in Section [2] We discuss the pipehne used to pro- 
cess the maps into an angular power spectrum in Sec- 
tion |3] We combine the SPT power spectrum with exter- 
nal data, most importantly the seven-year WMAP data 
release, to constrain cosmological models in Section |4] 
and we conclude in Section [5l 

2. OBSERVATIONS AND DATA REDUCTION 

The SPT is a 10-meter diameter off-axis Gregorian 
telescope located at the South Pole. The current receiver 
is equipped with 960 horn-coupled spiderweb bolome- 
ters with superconducting transition-edge sensors. The 
receiver included science-quality detectors at frequency 
bands centered at approximately 150 and 220 GHz in 
2008, and at 95, 150, and 220 GHz in 2009. The tele- 



scope and re c eiver are discussed in furt her detail in Ruhl 



pgn 



et al (20041, Padin et al. (2008), and Carlstrom et a. 



2.1. Fields and Observation Strategy 



spec- 
effec- 
The 



In this work we use data at 150 GHz taken dur- 
ing the 2008 and 2009 austral winters. This includes 
five fields whose locations, shapes, and effective areas 
(i.e. the area of the masks used in the power 
trum analysis) are given in Tabic [I] The total 
tive area is approximately 790 square degrees 
mean beam-convolved noise power in these fields is ap- 
proximated by the sum of a white noise component and 
a component that increases in power with decreasing £: 
d = (17.9 /iK-arcmin)2-|-3xlO-4(j^)-3 VK2|3The 
beam-deconvolved noise power and its two-component fit 
are shown in Figure [T] 

The fields were observed with two different types 
of scan strategies. The scan strategy used for the 
ra5h30dec-55 field consisted of constant-elevation 
scans across the field. After each scan back and forth 
in azimuth across the field, the telescope stepped 0.125° 
in elevation. We refer to a complete set of scans covering 
the entire field as an observation. 

The remaining four fields were observed using a 
lead/trail scan strategy. In this strategy each field was 
divided into two halves in Right Ascension, and the 
two halves were observed sequentially using constant- 
elevation scans. Due to the Earth's rotation, both halves 
of the field are observed at the same range of azimuth 
angle. This strategy allows for the possible removal of 
ground-synchronous signal, although we see no evidence 

The RA21HDEC-50 and RA21HDEC-60 fields overlap slightly. 
This reduces the effective total area of the power spectrum analysis. 
We have ignored this effect in our simulations, and have therefore 
underestimated the SPT bandpower errors by at most 0.4%. 

•^^ Throughout this work, the unit K refers to equivalent fluctu- 
ations in the CMB temperature, i.e., the temperature fluctuation 
of a 2.73 K blackbody that would be required to produce the same 
power fluctuation. The conversion factor is given by the derivative 
of the blackbody spectrum, evaluated at 2.73 K. 
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TABLE 1 
SPT Fields from 2008 and 2009 



Name 


R.A. (°) 


Decl. (°) 


AR.A. (°) 


ADecl. (°) 


Effective Area (sq. degrees) 


ra5h30dec-55 


82.5 


-55 


15.7 


10.3 


91.6 


ra23h30dec-55 


352.5 


-55 


18.2 


10.1 


105.5 


ra21hdec-60 


315.0 


-60 


30.5 


10.5 


156.9 


ra3h30dec-60 


52.5 


-60 


45.3 


10.6 


236.0 


ra21hdec-50 


315.0 


-50 


30.2 


10.5 


202.1 


Total 










792.1 



The locations and sizes of the fields observed by SPT in 2008 and 2009. For each field we give the center of the field in Right 
Ascension (R.A.) and Declination (Decl.), the extent of the field in Right Ascension and Declination, and the effective field 
area. 
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Fig. 1. — The beam-deconvolved noise power in the SPT maps 
used in this analysis (symbols show data and dotted line shows two- 
component model) compared to theoretical power spectra including 
CMB only [dashed line) and CMB+foregrounds [solid line). The 
precision of the power spectrum measurement is limited by sample 
variance rather than detector or atmospheric noise across most of 
the 650 < £ < 3000 range. 

for such a signal at the angular scales of interest and do 
not difference the lead and trail maps. 

2.2. Time-ordered Data to Maps 

Each SPT detector measures the sky brightness tem- 
perature plus noise, and records this measurement as 
the time-ordered data (TOD). The TOD are recorded 
at 100 Hz, so we have information in the TOD on sig- 
nals up to 50 Hz. For a typical scan speed and eleva- 
tion, 50 Hz corresponds to a mode oscillating along the 
scan direction a.t £ ^ 60, 000. This analysis, which only 
measures power at £ < 3000, can benefit computation- 
ally from using a down-sampled version of the TOD. We 
down-sample by a factor of six. Prior to down-sampling, 
we low-pass filter the TOD at 7.5 Hz. The combined 
effect of the filter and down-sampling is negligibly small 
(< 0.1% in power) on the scales of interest and we do 
not correct for it. 

The TOD are further low-pass filtered at 5 Hz as a 
safeguard against high-frequency noise being aliased into 
the signal band. The TOD are effectively high-pass fil- 
tered by the removal of a Legendre polynomial from the 
TOD of each detector on each scan. The order of the 
polynomial ranges from 9 to 18 and is chosen to have 



approximately the same number of degrees of freedom 
(dof) per unit angular distance (~1.5 dof per degree). 
The polynomial fit removes low-frequency instrumental 
and atmospheric noise. Regions of sky within 5 arcmin- 
utes of point sources with fiuxes S'isoghz > 50 mJy are 
masked during the polynomial fits. 

Correlated atmospheric noise remains in the TOD after 
the bandpass filtering. We remove the correlated noise 
by subtracting the mean signal across each bolometer 
wedgep^ at each time sample. This subtraction serves as 
an approximately isotropic high-pass filter. 

The data from each detector receives a weight based 
on the power spectral density of its calibrated TOD in 
the 1-3 Hz band, which corresponds approximately to the 
signal band of the power spectrum analysis presented in 
Section |3] We bin the data into map pixels based on 
the telescope pointing information. The maps us e the 
oblique L ambert equal-area azimuthal projection ( Sny- 



der 



1987) with 1' pixel resolution. The power spectrum 
ysis presented in Section [3] adopts the flat-sky ap- 
proximation, for which wavenumber k is equivalent to 
multipole £ and Fourier transforms replace spherical har- 
monic transforms. The 150 GHz map of the ra3h30dec- 
60 field is shown in Figure [2j 

2.3. Beam Functions 

The optical response as a function of angle — or beam — 
of the SPT must be measured accurately in order to cal- 
ibrate the signals in the maps as a function of angular 
scale. Due to the limited dynamic range of the detec- 
tors, the SPT beams were measured by combining maps 
of three types of sources: Jupiter, Venus, and the five 
brightest point sources in the CMB fields. The beam 
within a radius of 4' is measured on the five brightest 
point sources in these fields, and this naturally takes 
into account the enlargement of the effective beam due 
to random errors in the pointing reconstruction. Outside 
of this 4' radius, maps of Jupiter are used to constrain a 
diffuse, low-level sidelobe that accounts for roughly 15% 
of the total beam solid angle. Maps of Venus are used to 
stitch together the outer and inner beams. 

We measure the beam function, Bi, which is the az- 
imuthally averaged Fourier transform of the beam map. 
The beam function is normalized to unity at £ — 350. We 
express our uncertainty in the beam as an uncertainty in 
Bi. We account for the uncertainty arising from several 

'^'^ The SPT array consists of 6 wedge-shaped bolometer modules, 
each with 160 detectors. Wedges are configured with a set of filters 
that determine their observing frequency (e.g. 95, 150 or 220 GHz). 
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Fig. 2. — A map of the ra3h30dec-60 field, which is typical of the fields used in this analysis. The effective area is 236 square degrees. 
The structure visible in this map is due to primary CMB anisotropy, not instrumental or atmospheric noise. Modes with £ < 600 are 
strongly suppressed due to the hi gh-p ass filtering of the time-ordered data. The map has been multiplied by the apodization and point 
source masks described in Section |3.2[ such that bright point sources with SisoGHz > 50 mjy have been masked. A vertical stripe along 
the center of the map has been filtered more strongly than other regions. This stripe lies on the boundary of the lead and trail fields and is 
caused by high-pass filtering the time-ordered data by removing polynomial functions. This effect is accounted for in our analysis by using 
simulated observations. 
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statistical and systematic effects, including residual at- 
mospheric noise in the maps of Venus and Jupiter, and 
the weak dependence of Bg on the choice of radius used 
to stitch together the inner and outer beam maps. The 
different sources of uncertainty are incorporated into the 
power spectrum analysis through the bin-to-bin covari- 
ance matrix, as described in Section |3.6[ In Figure [3 
we show the beam functions and the quadrature sum of 
the different beam uncertainties, which gives a sense of 
the total uncertainty. There are small variations in the 
beam function between 2008 and 2009 across this multi- 
pole range due primarily to changes in the locations of 
the 150 GHz detector modules in the focal plane. The 
beam function is uncertain at the percent level across the 
multipoles of the power spectrum presented her e. The 
SPT b eams are discussed in more detail in .Lueker et al.l 
(2010) and Schaffer et al. (submitted). 
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Fig. 3. — The 150 GHz beam functions (bold, left scale) and frac- 
tional errors (thin, right scale). The beam function is normalized 
to one at £ = 350. 



2.4. Calibration 

The TOD are initially calibrated using a galactic HII 
region, RCW38. The final calibration used in this anal- 
ysis is calculated by comparing the SPT power spec- 
trum described in Section [3] with the seve n-year WMAP 
tempe rature power spectrum reported in Larson et al. 
(2011). The power spectra are compared across the an- 
gular range 650 < i < 1000, where the primary CMB 
anisotropy is the dominant source of power. First we 
construct a binned WMAP spectrum £)WMAP ^^^^^ jj^g^y 
be directly compared to the SPT spectrum by multi- 
plying the 6£ = 1 WMAP power spectrum by the SPT 
bandpower window functions described in Section |3.8[ 
Next we calculate the ratio R = Df^^^ /Df"^ aSd 
its weighted mean across this multipole range, (i?) = 



Y^i{wiRi) I Ylii'^i^ where Wi is the weight of the 



and is calculated as 



{D. 



WMAP\2 



bin 

>T ) . 



The WMAP error cr^jWMAP comes from the binned power 
spectrum provided by the WMAP team, which is also 
binned as Si = 50 at these multipoles. 
This calibration method is mo del- independent; it as- 



sumes only that the power measured in the SPT fields 
is statistically representative of the all-sky power mea- 
sured by WMAP. As such, this calibration method does 
not bias the constraints on cosmological parameters de- 
scribed in Section |4] We estimate the uncertainty in the 
SPT calibration to oe 3.1% in power. This uncertainty is 
included in the bin-to-bin covariance matrix as described 
in Section [XTI 



3. POWER SPECTRUM ANALYSIS 

In this section we describe the pipeline used to process 
the maps into an angular power spe ctrum. The method 
clos ely follows the approa ch used by Lueker et al. (2010) 



and ShirokofF et al 
imation, in which t 



(2011). We adopt the flat-sky approx- 
le angular wavenumber k is equivalent 
to multipole £ and Fourier transforms replace spherical 
harmonic transforms. The distortion to the power spec- 
trum due to adopting the flat-sky approximation on the 
SPT maps is negligibly small. We refer to the power in a 
given band of angular frequencies as the bandpower. We 
report bandpowers in terms of Di, where 



Df, = 



2n 



C, 



(1) 



3.1. Maps 

The power spectrum analysis begins with a set of 
150 GHz maps for each field. Each map corresponds 
to a single bottom-to-top observation of the field. For 
fields that were observed using a lead/trail method, we 
define a map to be a combination of consecutive lead and 
trail maps. We do not subtract the lead and trail maps. 
For the ra23h30dec-55 field, which was observed using 
large elevation steps, we take this one step further and 
define a map as the combination of a pair of consecutive 
lead/trail pairs. The composite map has more uniform 
coverage than the individual maps due to small elevation 
offsets between the individual maps. 

3.2. Windows 

The next step is to calculate the Fourier transform of 
each map, m"^, where A is the observation index. All 
maps of the same field are multiplied by the same win- 
dow W prior to the Fourier transform. This window 
is the product of an apodization mask used to avoid 
sharp edges at the map borders and a point source mask 
used to reduce the power from bright point sources. The 
apodization mask is a smoothed version of a map of ze- 
ros and ones in which a one denotes a pixel that was 
observed by at least one detector in every observation. 
We mask all point sources that we measure to have 
150 GHz flux > 50 mJy. Each point source is masked 
by a 5'-radius disk, with a Gaussian taper outside this 
radius with u taper = 5'. This relatively broad Gaus- 
sian taper was chosen to minimize the mode-coupling 
due to the point source mask. Using previous m easure 
ments of the mm-wave point sou rce population (Vieira 
et al. 2010 Shirokoff et al. 2011), we estimate that the 
power from residual point sources below this flux cut is 



This 



1.3 X 10-VK^ or ~ 18 ^JK^ {^] 
is approximately equal to the power from the primary 
CMB anisotropy at £ = 3000, the upper edge of the mul- 
tipole range of this analysis. A more aggressive point 
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source cut, say 10 mJy, could have been used to further 
reduce the residual power, but the gains were not consid- 
ered worth the cost of decreased sky area and increased 
mode-coupling. 

3.3. Cross-Spectra 

The next step is to cross-correlate maps from different 
observations of the same field. The noise in each individ- 
ual observation map is assumed to be uncorrelated with 
the noise in all other maps, so the resulting cross-spectra 
are free from noise bias. We calculate the cross-spectrum 
between maps from two different observations, A and B, 
and average within £-hms: 

5if^^/fc^F,i?e[m/mf*]\ , (2) 
\ 27r /^g^ 

where Hg is a two-dimensional weight array described 
below and is a vector in the two-dimensional, grid- 
ded Fourier plane of 5£ = 10 resolution. We average 
Df^ among all pairs of observations A and B, where 

A ^ B, to produce Di,. Each observation receives the 
same weight. For a typical field there are approximately 
200 maps and 20,000 cross-spectra. The maps are zero- 
padded prior to the Fourier transform such that the na- 
tive £ resolution is S£ = 10, which allows for clean sepa- 
ration into the final bins, which have a width of 6£ = 50. 
The lower edge of the lowest bin is £ = 650, while the 
upper edge of the upper bin is ^ = 3000. 

The noise in the maps used in this analysis is statisti- 
cally anisotropic; for a fixed £^ modes that oscillate per- 
pendicular to the scan direction (here defined as £x = 0) 
are the most noisy. For this reason, the modes that con- 
tribute to a given £-hin do not necessarily have uniform 
noise properties. We construct a two-dimensional weight 
array to optimally combine the modes contributing to 
each £-hm, 

He ^ {Cf- + N,)-\ (3) 

where Cf^ is the theoretical power spectrum used in the 
simulations described in Section [3 . 5 . 1 1 and Ni is the two- 
dimensional, calibrated, beam-deconvolved noise power. 
We smooth {Cf^ + Ng) with a Gaussian kernel of width 
(T£=425 in order to capture only the broad anisotropy of 
the noise power. The weight array is normalized such 
that J2eeb He. — \ for each bin h. The weights are ap- 
proximately uniform for the sample-variance-dominated 
bins {£ < 2700) and begin to preferentially de-weight 
£x ~ modes for bins at ^ > 2700. We estimate that 
this weighting scheme reduces our bandpower errors by 
approximately 8% in the highest bin. 

3.4. Jackknives 

Before proceeding with the rest of the power spectrum 
analysis, we apply a set of "jackknife" tests to the band- 
powers to search for possible systematic errors. In a 
jackknife test, the data are divided into two halves as- 
sociated with potential sources of systematic error. The 
two halves are differenced to remove any astronomical 
signal, and the resulting power spectrum is compared to 
an "expectation spectrum" , the power we expect to see 
in each jackknife spectrum in the absence of a systematic 



problem or spurious signal. This expectation spectrum 
can be non-zero due to mundane observational effects 
{e.g. a mid-season adjustment in the scanning strategy, 
or unequal weighting of left- and right- going scans). We 
use simulations to estimate the expectation spectra and 
find that they are small, with Di < 1 /iK^ at all mul- 
tipoles. Significant deviations of the jackknife spectrum 
from the expectation spectrum could indicate either a 
systematic contamination of the data or a misestimate 
of the noise. We construct difference maps (a single map 
in one jackknife half subtracted from a single map in 
the other half) and measure the jackknife spectrum as 
the average cross-spectrum between the difference maps, 
in a method analogous to that described in Section |3.3[ 
The jackknife spectra are calculated for five broad £-bins. 
We perform four jackknife tests based on the following 
criteria. 

• Time: We split the data into the first and second 
halves of observation. This tests for any system- 
atic that might be changing on weekly or monthly 
timescales. 

• Scan Direction: We split the data into left- 
going and right-going halves. This tests for scan- 
synchronous signals or any signal that is not time- 
symmetric, such as inaccurate deconvolution of the 
detector transfer functions. 

• Azimuthal Range: We split the data into observa- 
tions taken at azimuths that we expect to be more 
or less susceptible to ground pickup. We determine 
these azimuths by making maps using "ground- 
centered" coordinates (Azimuth/Elevation) as op- 
posed to "sky-centered" (R.A./Dec). Although we 
detect emission from the ground on large scales 
(£ 50) in these ground-centered maps, we do 
not expect such emission to bias our measurement 
of the sky power, as our observations are spread 
randomly in azimuth. We use the azimuth-based 
jackknife to test this assertion. 

• Moon: We split the data into observations taken 
at times when the Moon was either above or below 
the horizon. This tests for any significant coupling 
to the moon via far sidelobes of the SPT beam. 

For each test we calculate the of the jackknife spec- 
trum with respect to the expectation model. We cal- 
culate the probability to exceed (PTE) this for five 
degrees of freedom for each test and find PTE — 0.38, 
0.05, 0.92, 0.31 for the Time, Scan Direction, Azimuthal 
Range, and Moon tests, respectively. Although there 
are individual tests which have moderately low (0.05) or 
high (0.92) PTE's, they are consistent with a uniform 
distribution when taken as a whole. We therefore find 
no significant evidence for systematic contamination of 
the SPT bandpowers. 

3.5. Unbiased Spectra 

The spectra calculated in Section |3.3| are biased esti- 
mates of the true sky power. The unbiased spectra are 

Db^{K-'),,,D,, (4) 
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where b' is summed over. The K matrix accounts for the 
effects of the beams, TOD filtering, pixehzation, window- 
ing, and band-averaging. It can be expanded as 



mate. The factor W2 = J cfcW^ is a normaUzation factor 
for the area of the window. We then iterate on this esti- 
mate using the mode-couphng matrix: 



Kbb' = Pm (Mw [W] Fe, Bj, ) Qe,^' (5) 
where £ and £' are summed over. Qg and Pm are the 



binning and re-binning operators ( Hivon et al.||2002 ). 
Bi is the beam function described in Section |2.3[ J^'e s 
the transfer function due to TOD filt ering and map pix- 
ehzation, and is described in Section [3.5.11 The "mode- 
coupling matrix" Mu> [W] is due to observing a limited 
portion of the sky and is calculated analy tically f rom the 
known window W, as described in pueker et al. ( |2010] ). 
At the multipoles considered in this work, the elements 
of the mode-coupling kernel depend only on the distance 
from the diagonal. 



3.5.1. Transfer Function 

The transfer function Fi is calculated from simulated 
observations of 1500 sky realizations (300 per field) that 
have been smoothed by the appropriate beam. These 
simulations are also u sed to calculate the sample variance 
described in Section |3.6[ and it is therefore important 
that the power spectrum used to generate the simulated 
skies be consistent with previous measurements and with 
the power spectrum measured in this work. The simu- 
lated skies are Gaussian realizations of the primary CMB 
from the best-fit lensed WMAP-fCMB ACDM model 
from the seven-year WMAP releascp^ combined with con- 
tributions from randomly distributca point sources, clus- 
tered point sources, and the Sunyaev-Zel'dovich (SZ) ef- 
fect. The random point source component uses Ce = 
12.6 X 10-Vk' (i?^=30oo = Dsooo = 18.1 ^lK^). The 
clustered point source component uses Di — 3.5 jiK^fi, 
where fe = I ior i < 1500 and fi = for 
£ > 1500. This shape is designed to approximate the 
shape of the cl ustering power in both the linear and non - 
linear regimes fSh irokoff et"lI][20lT] [Millea et al.|[20T l| 



The SZ component uses the thermal SZ template ot |!je- 
hgal et al. (2010a I, which has a sha pe similar to the 
templates from more recent models ( Trac et al. 2011 



Shaw et al. 2010 Battaglia et al. 20lU), normalized to 
5.5 fiK at £ = 3000. The foreground compo- 



et al. 


(2011 


I and 


Vieira et al. 


(2010 



analysis. 

The simulated skies are observed using the SPT point- 
ing information, filtered identically to the real data, and 
processed into maps. The power spectrum of the sim- 
ulated maps is compared to the known inp ut spectrum 



to calculate the effective transfer function (Hivon et al. 



2002 1 using an iterative scheme. The initial estimate is 



F, 



(0) 



D, 



W2B/U 



,th ' 



(6) 



where the superscript (0) indicates that this is the first 
iteration in the transfer function estimates. We approxi- 
mate the coupling matrix as diagonal for this initial esti- 

■^^ http://lambda.gsfc.nasa.gov/ 



Mu'Flf^Be'^Df 



W2BtD 



2 nth 



(7) 



where £' is summed over. The transfer function estimate 
has converged after the second iteration, and we use the 
fifth iteration. The transfer function is equal to ~0.25 at 
£ = 650 and plateaus to ~0.85 at £ > 1200; the trans- 
fer function does not reach unity because of the strong 
filtering of £x ^ 300 modes. 

3.6. Bandpower Covariance Matrix 

The bandpower covariance matrix describes the bin-to- 
bin covariance of the unbiased spectrum and has signal 
term and a noise term. The signal term is estimated us- 
ing the bandpowe rs from the signal-only simulations de- 
scribed in Section [3. 5.1 and is referred to as the "sample 
variance." The noise term is estimated directly from the 
data using the distribution of the cross-spectrum band- 
po wers between o bservations A and B, as described 



Lueker et al. 



(2010), and is referred to as the "noise 
variance." TEe covariance is dominated by sample vari- 
ance at low multipoles and noise variance at high multi- 
poles, with the two being equal at ~ 2700. 

The initial estimate for the bandpower covariance ma- 
trix is poor for off-diagonal elements. We expect some 
statistical uncertainty. 



'^obs 



(8) 



This uncertainty is significantly higher than the true co- 
variance for almost all off-diagonal terms due to its de- 
pendence on the large diagonal covariances. We reduce 
the impact of this uncertainty by "conditioning" the co- 
variance matrix in the following manner. First we intro- 
duce the correlation matrix 



(9) 



The shape of the correlation matrix is determined by 
the mode-coupling matrix and is a function only of the 
distance from the diagonal. We calculate the conditioned 
correlation matrix by averaging the off-diagonal elements 
at a fixed separation from the diagonal: 



1 



(10) 



We set p'y = for all off-diagonal elements that are 
a distance t > 250 from the diagonal. The conditioned 
covariance matrix is then 



(11) 



We must also consider the bin-to-bin covariance due to 
the uncertainties in the beam function Bi as described in 
Section [2731 We construct a "beam correlation matrix" 
for each of the sources of beam uncertainty described in 
Section 1221 
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„6eam _ fSD,\ ( SDj 



where 



(12) 



(13) 



The combined beam correlation matrix is the sum of 
the beam correlation matrices due to each of the sources 
of uncertainty. The beam covariance matrix is then 



Jjeam 



(14) 



We calculate the beam covariance matrix for each year 
and add them to t he b andpower covariance matrix as 
described in Section [sItI 



3.7. Combining Different Fields 

We have five sets of bandpowers and covariances, one 
set per field, which must be combined into a single set of 
bandpowers and covariances. In the limit that the noise 
properties of all fields are identical, or in the limit that 
the precision of the power spectrum is limited by sam- 
ple variance on all scales of interest, each field would be 
weighted by its effective area (i. e. the area of its window) . 
While neither of these conditions is exactly true for our 
fields (the ra21hdec-50 field has higher noise than the 
other fields, and the power spectrum is dominated by 
noise variance at i > 2700), they are both nearly true, 
and we use the area-based weights. The field-averaged 
bandpowers and covariance are then 



where 



A' 



(15) 



(16) 



(17) 



is the area-based weight, and is the sum of the win- 
dow of the i"^ field. We introduce the beam covariances 
by first calculating the covariance for each year, then 
adding in the beam covariance for that year, and finally 
combining the covariances of the two years. The last step 
is to add the covariance due to the SPT calibration un- 
certainty, Clf — e^DiDj, where e = 0.031 corresponds 
to the 3.1% uncert ainty in the SPT power calibration 
discussed in Section [53 

The final bandpowers are listed in Table [2] and shown 
in Figure [4] 

3.8. Bandpower Window Functions 

In order to allow for a theoretical power spectrum Cf^ 
to be compared to the SPT bandpowers Cb, we calculate 
the bandpower window functions W^/l, defined as 



(18) 



Following the formalism described in Section |3.5[ we 
can write this as 



Cf" = {K-')hb'Pt'i'Me>iFeB^,Cf', 
which implies that 

W^,/£ = {K-\b'Pb't'M,.,F,Bl 



(19) 



(20) 



We calculate the bandpower window functions to be 
used with the final spectrum as the weighted average of 
the bandpower window functions from each field. 

4. COSMOLOGICAL CONSTRAINTS 

The SPT power spectrurrrj described in the previous 
section should be dominated by primary CMB anisotropy 
and can be used to refine estimates of cosmological model 
parameters. In this section we constrain cosmological 
parameters using the SPT power spectrum in conjunc- 
tion with data from the seven-ye ar WMAP data release 



(WMAP7, Larson et al. 



2oTTp 



and, 



m some cases, m 



conjunction with low-redshitt measurements of the Hub- 
ble co nstant i/o using the Hubble Space Telescope (Riessj 
et al.ir2011) and the baryon acoustic osc illation (BAOT 



featu re usmg SDSS and 2dFGRS data (Percival et al. 
20101 ). In the analyses that follow, the label "/^o+HAU" 
implies that the following Gaussian priors have been ap- 
plied: Ho = 73.8±2.4 km s'^ Mpc'^ rs/Dv{z = 0.2) = 
0.1905±0.0061; and rs/Dv{z = 0.35) = 0.1097±0.0036; 
where rg is the comoving sound horizon size at the baryon 
drag epoch, Dv{z) = [(1 + z)^D\{z)cz/H{z)]^^^ Da{z) 
is the angular diameter distance, and H{z) is the Hubble 
para meter. The inverse covariance matrix given in Eq. 



5 of jPercival et al. (|2010|) is used for the BAO measure- 



ments. 

4.1. Cosmological Model 

We fit the bandpowers to a model that includes four 
components: 

• Primary CMB. We use the standard, six- 
parameter, spatially flat, ACDM cosmological 
model to predict the power from primary CMB 
anisotropy. The six parameters are the baryon den- 
sity Vlbh? , the density of cold dark matter Clji^ , the 
optical depth of reionization r, the angular scale of 
the sound horizon at last scattering Bg, the ampli- 
tude of the primordial scalar fluctuations (at pivot 
scale fco — 0.002 Mpc^^) A|j,, and the spectral in- 
dex of the scalar fluctuations Us- The effects of 
gravitational lensing on the power spectrum of the 
CMB are calculat ed using a cosmology-de pendent 
lensing potential ( Lewis fc Challinor||2006 ). 



• "Poisson" point source power. Our model includes 
a term to account for the shot-noise fluctuation 
power from randomly distributed, emissive galax- 
ies. This term is constant in Ct and goes as 



Several of the data products presented in this work will be 
made available at http: / /pole.uchicago.edu/public/da ta/keislerll] 
and at http://lambda.gstc.nasa.gov/product/spt . 

•^^ We note that there is a small covariance between the SPT 
and WMAP bandpowers due to common sky coverage, but that it 
is negligibly small. The composite error has been underestimated 
by < 1% across the overlapping I range. 
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TABLE 2 

SPT Bandpowers and Bandpower Errors 



-^center 


De 


(j{Di) 


-^center 


Di 




■^center 


Di 


(j{Di) 


675 


1710 


95 


1475 


739 


25 


2275 


143 


5.4 


725 


2010 


98 


1525 


612 


21 


2325 


132 


5.1 


775 


2530 


110 


1575 


489 


17 


2375 


133 


5.0 


825 


2560 


110 


1625 


407 


14 


2425 


110 


4.5 


875 


2150 


93 


1675 


388 


14 


2475 


108 


4.4 


925 


1600 


69 


1725 


424 


13 


2525 


102 


4.7 


975 


1160 


51 


1775 


396 


13 


2575 


86.6 


3.8 


1025 


1100 


43 


1825 


343 


11 


2625 


83.6 


3.8 


1075 


1190 


46 


1875 


286 


9.8 


2675 


83.4 


4.1 


1125 


1250 


47 


1925 


234 


8.8 


2725 


83.8 


4.0 


1175 


1130 


43 


1975 


247 


8.2 


2775 


71.9 


3.6 


1225 


946 


35 


2025 


241 


8.2 


2825 


69.5 


3.9 


1275 


839 


29 


2075 


233 


7.7 


2875 


63.1 


4.1 


1325 


696 


27 


2125 


210 


7.0 


2925 


58.9 


3.5 


1375 


813 


27 


2175 


182 


6.2 


2975 


62.0 


4.1 


1425 


869 


28 


2225 


142 


5.4 









The SPT bandpowers and associated errors in units of iiK^. The 
errors do not include uncertainty in the SPT beam or cahbration. 



"Clustered" point source power. Our model in- 
cludes a term to account for the clustering of emis- 
sive galaxies. For this clustering contribution we 
use the template Df'" esc fg, where ft = I for 
£ < 1500, and ft = (y^)°-* for I > 1500. This 
shape is designed to approximate the shape of the 
clusterin g power in both the hnear and non-linea r 
regimes ( iShirokofr et al.||2011| [Millea et al.||2011[ ). 
We find that our cosmoiogicai constramts are not 
sensitive to the details of this shape. For example, 
there are no significant changes in the cosmoiogicai 
results if this shape is changed to a pure power law 



D 



CL 



SZ power. Our model includes a term to account 
for power from the thermal and kinetic SZ ef- 
fects. At the angular scales considered here, the 
two effects are expected to have similar shapes in 
^-space. We ther efore adopt the therm al SZ tem- 



plate provided in Sehgal et al. (2010a I, which has 



recent models ( 


Trac et al. 


2011 


Shaw et al. 


2010 


Battaglia et al. 


201U), to account for the total SZ 



power. 

For the purposes of this analysis, the primary CMB en- 
codes the cosmoiogicai information, while the last three 
components, the "foreground" terms, are nuisance pa- 
rameters. The foreground terms are used only when cal- 
culating the SPT likelihood; they are not used when cal- 
culating the WMAP likelihood. In our baseline model, 
we apply a Gaussian prior on the amplitude of each of 
the foreground terms. The prior on the Poisson power 
is -Dgg^QQ = 19.3 ± 3.5 /iK^ and is based on the power 
from sources w i th Si bqghz < 6.4 mJy, as measured in 
Shirokoff et al. (20111, and the power from sourc es with 
6.4 rt i Jy < Sib oghz < 50 niJy, as m easured in |Vieira| 
et al.| ( |2010[ ) and |Marriage et al.| ( |2011[) . The priors on the 
clustered power and SZ power are -D3Q00 = 5.0 ± 2.5 /xK^ 



and 



5.5 ± 3.0 /iK , as measured in Shirokoff 



et al. (2011). The widths of these priors span the mod- 
eling uncertainties in the relevant papers. Finally, we 
require the foreground terms to be positive. We find 
that our constraints on cosmoiogicai parameters do not 
depend strongly on these priors, as discussed below. 

One foreground that we have not explicitly accounted 
for is the emission from cirrus-like dust clouds in our 
Galaxy. Using a procedure similar to that described in 
Hall et al. (2010), we cross-correlate the SPT maps with 
predictfons tor the galactic dust emi ssion at 150 GHz 
in the SPT fields using model 8 of |Finkbeiner et ah] 
([1999 ). We detect the galactic dust in the cross- 
correlation and estimate the power from galactic dust 
in the field- averaged SPT spectrum to be — (1.4 ± 

0.4) (^)"^ Vk^. This is small compared to the SPT 
bandpower errors; subtracting this component changes 
the by 0.04 when all other parameters are kept fixed 
and by a negligible amount if the other foreground pa- 
rameters are allowed to move by amounts that are small 
compared to the widths of their priors. We conclude that 
galactic dust does not significantly contaminate the SPT 
power spectrum. 

The total number of parameters in our baseline model 
is nine: six for the primary CMB and three for fore- 
grounds. The nine-dimensional space is explored using 
a Markov Chain Monte Carlo (MCMC) technique. We 



usethe CosmoMCj^ software package ( Lewis & Bridle 
2OO2I , which itself uses the CAMg^ software package 



( [Lewis et al. 2000[) to calculate theTensed CMB power 
spectra, i'he CMi3 temperature and polarization spec- 
tra are calculated by CAMB for each cosmology. These 
spectra are passed to the likelihood software provided by 
the WMAP teanEi] to calculate the WMAP likehhood. 

http:/ /cosmologist.info/cosmomc/ 

32 http://camb.info/ . We use RECFAST Version 1.5. 

33 W e use Version 4.1 o f the WMAP UkcUhood code, avail- 
able at |http://lambda. gsfc.nasa.gov/ . In order to be consis- 
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Fig. 4. — The SPT power spectrum is shown in the left panel. The peak at ^ ~ 80 is the third acoustic peak. Fo r comparison we show 
in_the right panel other recent measurements of the CMB damping tail from ACBAR fReichard t et al.|2009[ l, QUaD JFriedman ct al.'2009j 

[Bro wn et al. 2009), ACT (Das et al. 2011a), and SPT (Shirokoff et al. 2011). The bandpower errors snowri in th ese p anels do not mclude 
beam or calibration u ncertainties. The ACT spectrum extends to ^ = 10,000. The previous SPT spectra, from [Lueker et al.| | |2010| ) and 

[Shirokoff et al.| | |2011| ), spanned the angular range 2000 < i < 10,000 and targeted secondary CMB anisotropy. 
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The SPT likelihood is calculated using the bandpowers, 
covariances, and bandpower window functions described 
in Section [Sj The age of the universe is required to be be- 
tween 10 and 20 Gyr, and the Hubble constant is required 
to be 0.4 < ft, < 1.0 where Hq = h 100 km s^^ Mpc"^ 
We assume that neutrinos are massless. 

Before combining the SPT and WMAP likelihoods, we 
first check that they independently give consistent con- 
straints on the six cosmological parameters. Because the 
scalar amplitude A|j and the optical depth r are com- 
pletely degenerate given only the high-i! SPT bandpow- 
ers, we fix r = 0.088 for the cosmological model that is 
constrained using only SPT data. We find that WMAP 
alone gives {lOOflbh^ .fi^h"^ , 10061s, n^, lO^Afj} = {2.24 ± 
0.056, 0.112 ± 0.0054, 1.039± 0.0027, 0.971 ± 0.014, 2.42 ± 
0.11}, while SPT alone gives {2.19 ± 0.18,0.110 ± 
0.013, 1.043 ± 0.0022, 0.953 ± 0.048, 2.49 ± 0.49}. Thus 
the two likelihoods are consistent with each other, and 
we proceed to combine them. 

The best-fit model from the joint SPT+WMAP like- 
lihood is shown in Figure [5] The baseline model pro- 
vides a good fit to the SPT data. The xVdof is 35.5/38 
(PTE=0.58) if the six cosmological parameters are con- 
sidered free and 35.5/44 (PTE=0.82) if the six cosmo- 
logical parameters are considered to be essentially fixed 
by the WMAP data. 

The marginalized likelihood distributions for the six 
cosmological parameters are shown in Figure |6] The ad- 
dition of the SPT data improves the constraints on f2b/i^ 
and rig by ~25% and the constraint on 6^ by nearly a 
factor of two. The parameter constraints for the base- 
line model are summarized in Table [3j and constraints 
on this model using SPT+WMAP-|-i7o+BAO are given 
in Table m 

The scalar spectral in dex n.s is le ss than one in sim- 
ple models of inflation (Linde 2008 ) . Recent measure- 
ments of Us come from WMAP7, = 0.967 ± 0.014 



dKomatsu et al.|201lD, ACBAR-FQ UaD-t-WMAP7, 
C).966;^:i;j3 dKomatsu et al.||201l|), and A CT-HWMAP7, 
n, = 0.962 ± 0.013 ( [Dunkley et al.| |2010[ ). The 

SPT+WMAP constraint IS 



0.9663 ±0.0112. 



(21) 



This is a 3. Oct preference for < 1 over the Harrison- 
Zel'dovich-Peebles index, ns= 1. This constraint is not 
significantly altered if we double the width of the priors 
on the foreground terms, in which case Us— 0.9666 ± 
0.0112. The constraint is also robust to doubling our 
uncertainties on the SPT beam functions or SPT cali- 
bration, in which cases ng— 0.9671 ± 0.0113 and Ug — 
0.9661 ± 0.0111, respectively. When the Hq and BAG 
data are included, the constraint is 71,= 0.9668 ± 0.0093, 
a 3.6cr preference for Us < I. 

4.2. Model Extensions 

In this section we consider extensions to our baseline 
model. These models continue to use a spatially flat, 
ACDM cosmological model, but allow a previously fixed 
parameter — the strength of gravitational lensing, the am- 
plitude of tensor fluctuations, the running of the spectral 

tent with our SPT-I-WMAP Markov chains, we recalculate all 
WMAP-only Markov chains rather than use those available at 
http://lambda.gsfc.nasa.gov/ . 



index, the primordial helium abundance, or the number 
of relativistic species — to vary freely. The structure of 
this section closely follows the clear presentation and dis- 
cussion of th e ACT+WMAP cons traints on parameter 
extensions by Dunkley et al. (2010), and therefore allows 
a straightforward comparison ot these similar datasets. 
We summarize constraints on these extension parame- 
ters using recent CMB datasets in Table [5j 

We note that we have also considered extensions with 
a free dark energy equation of state w or with massive 
neutrinos, and found that the addition of SPT data did 
not significantly improve upon the constraints on these 
models from WMAP alone. 

4.2.1. Gravitational Lensing 

The paths of CMB photons are distorted by the gravity 
of intervening matter as they travel from the surface of 
last scattering to us, a process referred to as gravitational 
lensing. The typical deflection angle is a few arcminutes, 
and the deflections are coherent over degree scales. Lens- 
ing encodes information on the distribution of matter at 
intermediate redshifts, and this information can be par 
tially recovered using correla tions of the GMB te mper 
ature and polarization fields ( Bernard cau 1997; Zaldar 
riaga fc Seljak|[T999l [Hu|[200l|]^^Gkamoto 2002| [ 



Previous etiorts to detect lensing ot the CJV^B include 
'^Sa detections from correlating quadratic reconstruc- 
tions of the lensing field from the GMB wit h other mass 
tracers ( [Smith et al.|2007||Hirata et al.|2008l ), < 3a detec- 
tions from the temperatiore pow^er spectrum (Reichardt 



et al.|2009[ jGalabrese et al.||2008[ [Pas et al.|2011ap , a 2a 



detection from the tour-pq int function ot the tempera- 
ture field (Smidt et al. 2011 ), and a 4a detecti on from the 
four-po int junction of the temperature field (Das et al. 
201 lb| . 

i'he GMB temperature power spectrum is altered by 
lensing at the few percent level. The acoustic peak 
structure is smoothed, and power is preferentially added 
to smaller angular scales. The SPT bandpowers accu- 
rately measure the acoustic peaks in the damping tail and 
should be sensitive to these effects. The effect of lensing 
on the GMB temperature power spectrum is largely cap- 
tured b y considering large-scale structure in the linear 
regime (Lewis & Ghallinor 2006). The corrections due 
to non- linear structures are small compared to the SPT 
bandpower errors and are not considered in this analysis. 

As a simple measure of the preference for lensing, we 
fit the SPT+WMAP bandpowers to a model in which 
the GMB is not lensed but that is otherwise identical 
to our baseline model. All parameters are free in this 
no-lensing model. We compare the likelihood of the 
best-fit no-lensing model to the likelihood of the best-fit 
lensed model. We find Ax^ = 2(ln/:i™'' - ln/:"° ^™'') = 
23.6, corresponding to a 4.9ct preference for the lensed 
model. This preference remains if we double the widths 
of the priors on the foreground parameters, in which case 
Ax^ — 23.0, or if we double the SPT beam uncertainty, 
in which case Ax^ = 22.9. We note that the distortion to 
the power spectrum caused by adopting the flat-sky ap- 
proximation on the SPT maps is small compared to the 
effect of gravitational lensing on the power spectrum. We 
estimate that correcting for this distortion would change 
the SPT likelihood by Ax^ < 0.04. 

To better quantify the strength of lensing preferred by 
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Fig. 6. — The one-dimensional marginalized constraints on the six cosmological parameters in the baseline model. The constraints from 
SPT+WMAP are shown by the blue solid lines, while the constraints from WMAP alone are shown by the orange dashed lines. 
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TABLE 3 

Constraints on Cosmological Parameters using SPT+WMAP 







A /~l T~\ TV T 


ACUM 


ACDM 


ACUM 


AODM 


ACUJVl 








A 

+ 


+r 


+ dus /d\n k 


+ Yp 


1 AT 


Primary 




2.22 ± 0.042 


2.22 ± 0.044 


2.24 ± 0.045 


2.21 ± 0.043 


2.26 ± 0.048 


2.27 ±0.054 
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n 1 1 o _i_ n nn/1 Q 
U.liz it U.UU4o 


U.liz it U.UUoU 


U.iuy it U.UUoU 


U.li / it u.uuoy 


U. 114 it U.UUOl 


n 1 OK -L- n ni o 

XJ.LZD It U.UIZ 






1.041 ± 0.0016 


1.041 ± 0.0016 


1.041 ± 0.0016 


1.041 ±0.0017 


1.043 ±0.0020 


1.040 ± 0.0019 




ris 


0.9663 ±0.0112 


0.9655 ±0.0114 


0.9743 ±0.0128 


0.9732 ±0.0120 


0.9793 ± 0.0140 


0.9874 ± 0.0193 






0.0851 ± 0.014 


0.0853 ± 0.014 


0.0860 ± 0.014 


0.0909 ± 0.016 


0.0884 ± 0.015 


0.0883 ±0.015 






2.43 ±0.10 


2.44 ± 0.11 


2.36 ±0.11 


2.37 ± 0.11 


2.39 ±0.10 


2.37 ± 0.11 


Extension 


^0.65 




0.94 ± 0.15 










Parameters 


r 






< 0.21 










dus/dlnk 








-0.024 ±0.013 








Yp 


(0.2478 ± 0.0002) 


(0.2478 ± 0.0002) 


(0.2479 ± 0.0002) 


(0.2477 ± 0.0002) 


0.296 ±0.030 


(0.2579 ± 0.008) 




N,ti 


(3.046) 


(3.046) 


(3.046) 


(3.046) 


(3.046) 


3.85 ±0.62 


Derived 




(0.814 ±0.024) 


(0.816 ±0.025) 


(0.805 ±0.025) 


(0.832 ±0.027) 


(0.837 ±0.029) 


(0.859 ± 0.043) 






7506.5 


7506.2 


7506.4 


7503.6 


7504.4 


7505.5 



The constraints on cosmological parameters using SPT+WMAP7. We report the mean of the Ukelihood distribution and the symmetric 68% 
confidence interval about the mean. We report the 95% upper limit on the tensor-to-scalar ratio r. Parameters labeled with " — " were held 
at their default values: Aj^ = l,r = 0, and dus/dlnk = 0. 



TABLE 4 

Constraints on Cosmological Parameters using SPT+WMAP+Hq+BAO 







ACDM 


ACDM 


ACDM 


ACDM 


ACDM 


ACDM 








+ Al 


+r 


+ dris/dln k 


+ Yp 


+ iVcff 


Primary 


lOO^j/i^ 


2.23 ±0.038 


2.22 ±0.039 


2.24 ± 0.040 


2.23 ±0.040 


2.27 ± 0.044 


2.26 ±0.042 


Parameters 




0.112 ±0.0028 


0.112 ±0.0029 


0.112 ±0.0030 


0.114 ± 0.0031 


0.114 ±0.0032 


0.129 ±0.0093 




1009s 


1.041 ± 0.0015 


1.041 ±0.0016 


1.041 ± 0.0015 


1.041 ± 0.0016 


1.043 ±0.0020 


1.039 ±0.0017 




ris 


0.9668 ± 0.0093 


0.9659 ± 0.0095 


0.9711 ± 0.0099 


0.9758 ±0.0111 


0.9814 ±0.0126 


0.9836 ±0.0124 




T 


0.0851 ± 0.014 


0.0852 ±0.014 


0.0842 ± 0.014 


0.0934 ± 0.016 


0.0890 ± 0.015 


0.0859 ±0.014 






2.43 ±0.082 


2.44 ± 0.085 


2.39 ±0.088 


2.35 ±0.095 


2.39 ±0.085 


2.41 ± 0.084 


Extension 






0.95 ±0.15 










Parameters 


r 






< 0.17 










dris/dlnk 








-0.020 ±0.012 








Yp 


(0.2478 ± 0.0002) 


(0.2478 ± 0.0002) 


(0.2478 ± 0.0002) 


(0.2478 ±0.0002) 


0.300 ±0.030 


(0.2581 ± 0.005) 






(3.046) 


(3.046) 


(3.046) 


(3.046) 


(3.046) 


3.86 ±0.42 


Derived 




(0.818 ±0.019) 


(0.818 ±0.019) 


(0.816 ±0.019) 


(0.824 ± 0.020) 


(0.841 ± 0.024) 


(0.871 ± 0.033) 






7510.7 


7510.6 


7510.7 


7507.8 


7508.0 


7507.4 



The constraints on cosmological parameters using SPT+WMAP7+i?o+BAO. We report the mean of the likelihood distribution and the 
symmetric 68% confidence interval about the mean. We report the 95% upper limit on the tensor-to-scalar ratio r. Parameters labeled with 
" — " were held at their default values: Aj^ = l,r = 0, and dris/dlnk = 0. 



TABLE 5 

Constraints on Model Extensions using Recent CMB Datasets 





WMAP7 


ACBAR+QUaD+WMAP7 


ACT+WMAP7 


SPT+WMAP7 


r 


< 0.7 


< 0.33 


< 0.25 


< 0.21 


dus/d In k 


[-0.084, 0.020] 


[-0.084, 0.003] 


-0.034 ± 0.018 


-0.024 ± 0.013 


Yp 


< 0.51 


0.326 ± 0.075 


0.313 ±0.044 


0.296 ±0.030 




> 2.7 




5.3 ± 1.3 


3.85 ± 0.62 



The constraints on cosmological parameters in certain model extensions using recent CMB 
data sets. We use WMAP7 (Larson et al. 2011, Komatsu et al. 2011), ACBAR (Reichardtl 
I et al.|2009| ), QUaD |,Brown et al.|2009| ), ACT I jDas et al.|2011a| ), and SPT (this work). AIT 
upper and lower limits and all two-sided limits [shown m brackets) are 95%. 
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the data, we consider the parameter Al which rescales 
the lensing potential power spectrum, 



(22) 



As in all of our models, Cf''' is calculated in a cosmology- 
dependent manner at each point in the Markov chain. 
The Ai^ parameter is phenomenological and we allow it 
to be negative. In the standard scenario, Al = 1, while 
Al = corresponds to no lensing. 

This parameter has been constrained using measure- 
ments of the CMB damping tail in conjunction with 
WMAP data, and more recently by measuring the lens- 
ing signal encoded in the CMB te mperature four-point 
function. iReichardt et all (|2009|) used A CBAR data 




in conjunction with hve-year WMAP d ata (Nolta et al 

while Calabrese et al. 
found = 3.0±0.9 using the same datasets. Das 
(|2011a|) used power spect ra from ACT-HWM AP7 



to measure A^ ~ 1.3 ± 0.5, and Das et al. (2011b) used 
the four-point function of the AC'l' temperature maps to 
measure Al = 1.16 ±0.29. 

The constraints on A^ from SPT-I-WMAP7 are shown 
in Figure [7j While the constraint on Al is non-Gaussian 
in shape, we find that the constraint on A^^^ (which still 
has an expectation value of 1) is approximately Gaussian. 
With SPT-I-WMAP7, we finc^ 



A 



0.65 



= 0.94 ±0.15. 



(23) 



This constraint does not change significantly if we double 
the width of the priors on the foreground terms, in which 
case = 0.94 ±0.15. The SPT+WMAP7 data reject 

a non-lensed CMB and are consistent with the expected 
level of lensing. This provides a consistency check on 
the standard picture of large-scale structure formation. 
Ongoing and future measurements of CMB lensing will 
move beyond this consistency check and constrain pa- 
rameters that affect the growth of large-scale structure, 
such as the prop erties of dark energy and the sum of th e 
neutrino masses ( Smith et al.|2006 Sherwin et al.|20lT ). 



4.2.2. Tensor Perturbations 

Infiation is expected to produce primordial tensor per- 
turbations (i.e. gravitational waves). These perturba- 
tions imprint potentially detectable effects onto the CMB 
temperature and polarization spectra. The amplitude of 
the tensor spectrum is often given in terms of the the 
tensor-to-scalar ratio, r = A^(/co)/A^(A:o) with pivot 
scale fco = 0.002 Mpc~^j^ A detection of r would pro- 
vide an extremely interesting window onto the energy 
scale of inflation. 

Measurements of the B-mode polarization at low multi- 
poles will ultimately provide the strongest constraints on 
r. To date, the best constraint on r from B-mode polar- 
ization comes from the BICEP experiment ( Chiang et al. 



Note that this constraint, along with the other constraints 
on A^^^ listed in this work, impHcitly assumes a uniform prior 
on Al. However, the result does not change significantly if we 
modify the prior to be uniform in A^^^ instead, in which case 

A°f^ =0.93 ±0.15. 

We assume that the spectral index of the tensor perturbations 
is nt = — r/8. 



20101, giving r < 0.7 (95% CL). Stronger constraints 



are currently placed using WMAP's me asurement of the 
temp erature and polarization spectra (Komatsu et al. 
2011] ), which give r < 0.36 (95% CL). 

The constraint on r can be improved indirectly with 
the addition of small-scale CMB measurements. The 
CMB power at low multipoles increases as r increases, 
but this effect can be partially cancelled by increasing 
Us and decreasing A^. The small-scale CMB measure- 
m ents help to break these degeneracies, as demonstrated 



Komatsu et~a l] (|2011|, which found r < .33 (95% 



CL) u s ing A rnAR+QUaD+WMAP7, and in Dunkley 



et al. (2010), which found 
XCT±WMAP7. The SPT 
to be 



r < 0.25 (95% CL) using 
fWMAP7 data constrain r 



< 0.21 (95% CL). 



(24) 



When the Hq and BAO data are added, the constraint 
improves to r < 0.17 (95% CL). Figure Is] shows the one- 
dimensional marginalized constraint on r and the two- 
dimensional constraint for r and the spectral index rig. 
We show the predi ctions for r and from chaotic in- 
flationary models (Linde 1983) with inflaton potential 
V{(f)) oc (pP and N = 60, where TV is the number of e- 
folds between the epoch when modes that are measured 
by SPT and WMAP exited the horizon during inflation 
and the end of inflation. These models predict r = Ap/N 
and Ug — 1 — ijj + 2)/2N . Models with p > 3 are disfa- 
vored at more than 95% confidence for N < 60. 

4.2.3. Running of the Spectral Index 

The power spectrum of primordial scalar fluctuations 
is typically parametrized as a power law, 

k " 



AUk) = Afs(fco) 



fcn 



(25) 



In this section we allow the power spectrum to depart 
from a pure power law. The "running" of the spectral 
index is parametrized atPj 



A|(fc) = A|(fco) 



fco 



Tis — 1+ ^ \n(k/ko)d7is /d\nk 



(26) 



The running parameter dns/dlnk is predicted to be 
small by most inflationary theories, and a detection of 
a non-zero dns/dlnk could provide information o n the 
inflationary potential (Kosowsky & Turner 1995). Re- 
cent CMB constraints on the running include —0.084 < 



dnjdlnk < 0.020 (95% CL) from WMAP7 (Komatsu 
et al.||2011|), -0.084 < dnjdln k < 0.003 (95% CL ) 



from ACBAR+QUaD+WMAP7 ( Komatsu et al1|20TT ) 
and dnjdlnk = -0.034 ± 0.018 from ACi'+WlVlAPV 



(Dunkley et al. 2010). The SPT+WMAP7 data con- 



strain dns/dlnk to be 



dns/d In k = -0.024 ± 0.013. 



(27) 



The data mildly prefer, at 1.8cr, a negative spectral run- 
ning. The constraint is dns/dlnk = -0.020 ± 0.012, 

The factor of 1/2 is due to considering the effective change in 
ris in din A|j/c( In fc. 
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Fig. 7. — The one-dimensional marginalized constraint on the gravitational lensing parameter A^. This parameter rescales the lensing 
potential power spectrum as cf''' — > A^cf'^. 




Fig. 8. — The one-dimensional marginalized constraint on the tcnsor-to-scalar ratio r (left) and the two-dimensional constraint on r 
and the spectral index (right). The dashed line shows predictions for r and from chaotic inflationary models with inflaton potential 
V{4>) oc (pP and = 60, where A'' is the number of e-folds between the epoch when modes that are measured by SPT and WMAP exited 
the horizon during inflation and the end of inflation. The two-dimensional contours show the 68% and 95% confldence intervals. 
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a 1.7a preference for negative running, when the Hq 
and BAO data are added. As discussed in Section 14.3 
and shown in Table |6j the constraint is dus/dlnk — 
— 0.017 ± 0.012, a 1.4(7 preference for negative running, 
when information from local galaxy clusters is added. 
Figure [9] shows the one-dimensional marginalized con- 
straint on dus/dhik and the two-dimensional constraint 
for dng/dlnk and the spectral index Tigp^ 

4.2.4. Primordial Helium Abundance 

When the universe cools to T ~ 0.1 MeV, light nuclei 
begin to form , a process known a.s big bang nucleosyn - 
thesis (BBN) ( [Schramm fc Turner|1998||Steigman|2007l ). 
The primordial abundance (mass fraction) of ''He is de- 
noted as Yp and is a functio n of baryon density and t he 
expansion rate during BBN ( Simha & Steigman l[2008t : 



yp = 0.2485-F0.0016[(273.9fib/i2-6)-f 100(5-1)], (28) 



where 



= 1 + (7/43)(A^off - 3.046). 



(29) 



The factor generically accounts for any non- 

standard expansion rate prior to and during BBN, here 
parametrized in terms of the effective number of rela- 
tivistic particle species Ncs. We calculate Yp in this 
"BBN-consistent" manner in all of our models, unless 
noted otherwise!^ 
For our baseline model, we find Yp = 0.2478 ± 0.0002. 

[391 



This very small uncertaintjj^ on Yp is due only to the 
uncertainty in Q,bh?, as A'efr is fixed to its standard value 
of 3.046. This constraint relies on standard BBN theory 
being correct. We can consider an independent mea- 
surement of Yp-. the value of Yp preferred by the CMB 
due solely to the effect of helium on the CMB damp- 
ing tail. Helium combines earlier than hydrogen, and 
thus more helium (at fixed baryon density) leads to fewer 
free electrons during hydrogen recombination. This, in 
turn, leads to larger diffusion lengths for photons and 
less power in the CMB damping tail. 

A simple test of the preference of the CMB data 
for non-zero primordial helium follows. We compare 
the maximum likelihood in a model with no ^He to 
the maximum likelihood in our baseline model. Us- 
ing SPT-hWMAP7, we find that the standard, BBN- 
consistent helium abundance is preferred over no helium 
at 7.7cr (Ax^ = 58.8). 

The estimates for dus/dlnk and ng are highly corr elated for 
the ty pical pivot scale, fco = 0.002 Mpc~^. As in Dunk ley et al.| 
l |2010[ l, we have calculated the spectral index at a new, less corre- 
lated pivot scale kg = 0.015 Mpc~', where na(ko = .015 Mpc~^) 
= nJko = 0.002 Mpc-i)-|-ln(0.015/0.002)dns/dlnA: j Cortes et at] 
[2007| |. 

We note that Yp is calculated within CAMB using the 
PArthENoPE BBN code (Pisanti et al. 2008), and that the Yo 
calculated with this code ditters trom the tunction given in Eg. |28| 
by < 1% in the range Clbh^ = [0.021,0.023] and N^g = [2,6]. We 
pass ATVj. = (TVeff - 3.046) to PArthENoPE, such that our def- 
inition of A'efr reisers to N^ff in the epoch after electron-positron 
annihilation at T ~ 0.5 MeV. 

The ~0.1% uncertainty quoted here is the statistical uncer- 
tainty and is slightly smaller than the 0.2% theoretic al uncertainty 
on Y„ in the PArthENoPE code used in CAMB jPisanti et al.| 
[2008| . 



We can extend this test by promoting Yp to a free 
par ame ter. In such a model, Yp no longer obeys 



Eg. 28 but rather is free to vary. 
(|20iyusc WMAP7 to infer Yp < 0.51 



Komatsu et 
95% CL) 



al 
and 



ACH AR+QUaD-hWMAP7 to measure Yp = 0.326 ± 



0.075. 

sure 
strain Yr, to be 



Dunkley et al. (2010) use ACT 

0.313 ± 0.6k: 



WMAP7 to mea- 
'he SPT+WMAP7 data con- 



Fp = 0.296 ± 0.030. 



(30) 

The data mildly prefer, at l.Gtr, a value of Yp that is 
larger than the value we obtain assuming standard BBN 
theory. The constraint is Yp = 0.300 ±0.030, 1.7cr higher 
than the standard BBN value, when the Hp and BAO 
data are added. As discussed in Section 14.31 and shown 
in Table [6| the constraint is Yp = 0.288 ± 0.029, 1.4cr 
higher than the standard BBN value, whe n in formation 
from local galaxy clusters is added. Figure 10 shows the 
one-dimensional marginalized constraint on Yp. 

The primordial *He abundance may also be inferred 
from observations of low-metallicity extragalactic HII re 
gions (Izqtov e t al. 2007; Pcimbcrt ct al. 20Wj [Izotov 
fc Thu ahlUTU^ jAver et al.||2010n20n| . For example. 
Aver et al. (2011 ) provides an extensive analysis of the 
systematic uncertainties associated with these measure- 
ments and finds Yp = 0.2609 ±0.0117 (or 0.2573l:o°ooll if 
the metallicity slope dYp/dZ is required to be positive). 
These values lie between, and are consistent with, the re- 
suh from our baseline model, Yp = 0.2478 ± 0.0002, and 
the result from our free-Yp model, Yp — 0.296 ± 0.030. 



4.2.5. Number of Relativistic Species 

In the standard theory of the early universe, there are 
three neutrino species that contribute ~10% of the en- 
ergy density at recombination. The effective number of 
particle species that are relativistic prior to and during 
recombination, A"offj^ is slightly higher (3.046) due to 
energy injection from electron -positron annihila t ion at 
the end of neutrino freeze-out | D icus et al.j 1982 [Lopez 
et al. [1999 Mangano et al. 200S ). A signincant detec- 
tion of iVoff ^ 3.046 could point to the presence of extra 
relativistic species in the early universe. 

The addition of extra relativistic species increases the 
expansion rate during the radiation-dominated era. If 
the parameters that are robustly measured by WMAP — 
Qbh^, zeq and 0s — are held fixed, then the main effect 
of this increased expansion rate is to increase the angu- 
lar scale of photo n diffusion and thereby lower power in 



the damping tail (Hu fc White 1996 Bashinsky fc Sel ^ 
jak[[2004 Hou et al.|201ip . Conversely, measurements of 



the CMB damping tail can, in conjunction with WMAP, 
constrain the number of relativistic species. 

A simple test of the preference of the CMB data for a 
non-zero number of relativistic species follows. We com- 
pare the maximum likelihood in a model with Ncs=0 
to the maximum likelihood in our baseline model. Us- 
ing SPT+WMAP7, we find that the standard value of 

■^O TVcff is defined such that = Af<,ff7/8(4/ll)"'/3/3.^ 
*' ^EQ is the redshift at which matter and radiation have equal 
energy densities and satisfies (l+ZEQ)=f2m/f2r, where Qm is the 
density of matter and Qr is the density of radiation (photons and 
neutrinos). 



17 



1 .0 


//V 








SPT+WMAP 


0.8 


' \^ 


WMAP 




' \ ^ 




8 0.6 


/ / \ 


\ 




/ / \ 
/ / \ 


\ 


0.4 


' 1 \ 


\ 

\ 






\ 


0.2 


t 1 ' 


. \ 






\ \ 


0.0 


; / / 




-0 


10 -0.05 0.00 0.05 

dn^/ dink 



o 



LO 

q 
d 



.00 



0.95 



SPT+WMAP 
WMAP 




0.90 

-0.10 -0.05 0.00 0.05 

drig/ dink 



Fig. 9. — The one-dimensional marginalized constraint on the running of the spectral index dUg / d\r\k (left) and the two-dimensional 
constraint on dus / d\i\k and the spectral index rts (right). The two-dimensional contours show the 68% and 95% confidence intervals. 
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A'eff =3.046 is preferred over no relativistic species at 7.5a 
(Ax^ = 56.3). The CMB data strongly prefer the exis- 
tence of neutrinos over no neutrinos. 
We can ex tend this test by promo ting A^cff to a free 



parameter. Komatsu et al.| |2011') find Npff > 2.7 



fi% CL) using WMAP7 alone, while Dunkley et al. 
)10) find A^eff = 5.3 ± 1.3 using ACT+WMAP7. The 
T+WMAP7 data constrain N^g to be 



A^off = 3.85 ± 0.62. 



(31) 



This constraint is 1.3tT higher than the standard Aog = 
3.046. When the Hq and BAO data are added, the con- 
straint improves to A'cfi = 3.86 ± 0.42, 1.9(7 higher than 
the standard N^s = 3.046. As discussed in Section [473| 
and shown in Table[6) the constraint is A^eff = 3.42±0.32, 
1.2(7 higher than the standard value, when information 
from local galaxy clusters is added. Figure [TT] shows the 
one-dimensional marginalized constraint on A'off. 

Acff has also been constrained using measurements of 
abundances of ^He and deuterium. These abundances 
are sensitive to the expansion rate during BEN, which, 
in turn, is se nsitive to N^s. Using these methods, |Simha| 
fcSteigman| ( |2008[ ) find N^s = 2.4±0.4, ahhough the^He 
abundance used in that work, Yp — 0.240±0.006, is lower 
than more r ecent determinatio ns (e.g. Yp = 0.2609 ± 
0.0117 from |Aver et al] ( |2011| ). Mangano fc Serpic0| 
(2011) use *He and deuterium abundances to provide a 



conservative upper limit of A^cff< 4.2 (95% CL). 

The effective number of relativistic species has also 
been constrained by combining low-redshi ft measure- 
ments with WMAP5 data. For example, |Reid et al.| 

using the abundance of 



(2010) found A^off = 3.76 



-1-0.63 
-0.68 



optically-selected galaxy cluste rs, CMB data, and a mea- 
surement of Hq. Similarly, iMantz et al. (2010a I found 



A"off — 3.4_g'5 using the abundance of X-ray-selected 
galaxy clusters, galaxy cluster gas mass fraction data, 
CMB data, supernova data, BAO data, and a mea- 
surement of Hq. These results are consistent with, but 
not independent from, the constraints presented here, as 
WMAP data is common to all of these constraints. 

4.3. Discussion of Models that allow for Additional 
Damping 

In the previous sections we have fit the SPT-f WMAP 
data (and in some cases, Hq and BAO data) to a flat, 
ACDM cosmological model and to a number of exten- 
sions beyond this model. Of these extensions, three of 
them — allowing for spectral running, varying the pri- 
mordial helium abundance, and varying the effective 
number of relativistic species — improved the fit to the 
data by Ax^ ~ 3, primarily by lowering the predicted 
power in the CMB damping tail. In other words, mod- 
els with either a negative spectral running, a high he- 
lium abundance, or a high effective number of relativistic 
species were mildly preferred, at 1.6-1.9ct depending on 
the model. We use this section to briefly discuss the ro- 
bustness of the preference, the degeneracies between the 
parameters in these models, and the consistency (or lack 
thereof) of these models with external datasets. We find 
that these models prefer values of as that are disfavored 
by measurements of local galaxy clusters, and that in- 
cluding the cluster information brings the constraints on 
these parameters closer to their standard values. 



4.3.1. Robustness of Preference 

The models with spectral running and free Yp prefer 
non-standard models at l.Scr and 1.6cr, respectively, us- 
ing only the CMB data, and these preferences shift to 
1.7(7 if the Hq and BAO data are included. The model 
with free A^off deviates from the standard Aoff = 3.046 
by 1.3(7 when only CMB data are used, and a high Ncs 
is preferred at 1.9(7 only after the Hq and BAO data are 
combined with the CMB data. The data do not show any 
significant preference for one extension over the others. 

We have tried doubling the SPT beam uncertainty, 
doubling the widths of the priors on the foreground pa- 
rameters, and using a ''Df^ = constant" template for 
the clustered point source power, and find that none of 
these significantly weakens the preference for additional 
damping. 

4.3.2. Degeneracies between Parameters 

The best-fit model spectrum is lower by ~2.5% at 
£ = 2000 for each of the three models relative to the 
baseline model. This suggests that the extension param- 
eters dus/dlnk, Yp, and A^eff are degenerate. Indeed, 
we find that for models in which two or three of these 
parameters are free, there are degeneracies between the 
extension parameters, and the constraint on any one pa- 
rameter is weakened. 

To illustrate this degeneracy, consider a model in which 
the primordial helium abundance Yp and the effective 
number of relativistic species N^s are both free. In 
this model, the two are no longer related by BBN the- 
ory and are independent. In Figure |12[ we show the 
two-dimensional constraint on these parameters using 
SPT-I-WMAP7. First, note that there is a degeneracy 
between the parameters; a cosmology with a high A^eff 
can be accommodated by lowering Yp. But the two are 
not completely degenerate, as the contours do not ex- 
tend to Acff^O or Yp—0. The marginalized constraints 
are A^eff = 3.4 ± 1.0 and Yp = 0.283 ± 0.045. This is an 
interesting result in itself: the SPT-I-WMAP7 data are 
able to significantly detect the effects of helium and neu- 
trinos independently. The main point, however, is that 
the three extension parameters that affect damping-tail 
power — dus/d In k,Yp, and Nes — are degenerate, and the 
constraint on any one of them is reduced if we allow the 
others to be free. For simplicity we have presented the 
results from models where only one of these parameters 
is free, rather than from the models where two or more 
of them are simultaneously free. 

4.3.3. Consistency with External Data 

Are any of these scenarios — a negative spectral run- 
ning, a high Yp, or a high A'cft — obviously ruled out 
by existing data? The spectral running is primarily 
constrained through the CMB, and previous measure- 
ments do not rule out the SPT-HWMAP7 best-fit value, 
dus/dlnk = -0.024 ± 0.013. Regarding helium, the 
SPT+WMAP7 best-fit value of Yp = 0.296 ± 0.030 is 
consistent at 1.2a wi t h the Hll-region-based measure 



ments of Aver et al. (2011), who find K = 0.2609 ± 



0.0117. Regardmg the number of relativistic species, the 
SPT+WMAP7-hi/o+BAO best-fit value for the effective 
number of relativistic species, A"eff = 3.86 ± 0.42, is con- 
sistent with the upper hmit of A^eff < 4.2 (95% CL) 
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Fig. 11. — The one-dimensional marginalized constraint on the effective number of relativistic species N^g. The standard value of 
N^ff = 3.046 is shown by the vertical dotted line. 



TABLE 6 

Constraints on Cosmological Parameters using 
SPT+WMAP+Ho+BAO+Clusters 







ACDM 


ACDM 


ACDM 






+ dus/dln k 


+ Yp 


+ A'off 


Primary 


lOOQft/i^ 


2.23 ±0.040 


2.26 ±0.045 


2.24 ±0.041 


Parameters 




0.111 ± 0.0020 


0.111 ± 0.0020 


0.116 ±0.0054 




1006^ 


1.041 ± 0.0016 


1.043 ±0.0019 


1.040 ±0.0017 




rig 


0.9751 ± 0.0110 


0.9787 ± 0.0123 


0.9757 ± 0.0116 






0.0897 ± 0.015 


0.0852 ± 0.014 


0.0821 ± 0.014 




109A2„ 


2.33 ±0.092 


2.35 ±0.082 


2.37 ±0.081 


Extension 


dus/dlnk 


-0.017 ± 0.012 






Parameters 


Yp 


(0.2478 ± 0.0002) 


0.288 ±0.029 


(0.2526 ±0.004) 




N,a 


(3.046) 


(3.046) 


3.42 ±0.32 


Derived 


OS 


(0.809 ±0.014) 


(0.819 ±0.016) 


(0.823 ±0.019) 




■ 


7509.3 


7509.3 


7510.3 



The constraints on cosmological parameters using 

SPT+WMAP7+//o+BAO+Clusters, where "Clusters" refers to the local 
cluster abundance measurement of Vikhlinin et al. (2009). We report the mean 
of the likelihood distribution and the symmetric 68% confidence interval about 
the mean. The label " — " signifies dus/dlnk = 0. 
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inferred by Mangano & Serpico (2011) using ''He and 
deuterium abundances. It also c onsistent wi t h, bu t not 
independent from, t he results of Reid et aT] (2010) and 
Mantz et aL] ( |2010a| ), which combine WMAF5 data with 
low-redshift measurements. 

A common feature of these three models is that they 
require values of as, the amplitude of linear matter 
fluctuations on scales of 8 Mpc a,t z — 0, that 
are higher than those favored in the baseline model. 
The as required in the high-A'eff model is particularly 
high: the constraint on erg is as = 0.871 ± 0.033 us- 
ing SPT-hWMAP+i/o+BAO, while the equivalent con- 
straint for the baseline model is as = 0.818±0.019. This 
correlati on between high-A^eff models and high erg was 
noted in Dunkley et al. ( 2010[ ). Are such high values of 
as consistent with low-redshift measurements? To an- 
swer this question, we consider the g alaxy cluster abun- 
dance measurement of Vikhlinin et al.^ ( ,2009) , which di- 
rectly and tightly constrains erg and is consistent with 
other measurements of structure at low to medium red- 



shift dRozo et al.||2010 , Mantz ct al.| 


2010bl 


Vanderlinde 


et al. 201O| |Sehgal et al.i 


2010bp. The aut 


tiors use the 



fer CT8(r2™/0.25)0-4^ = 0.813 ± 0.013 ± 0.024, where the 
second set of errors is an estimate of the systematic un- 
certainty due to the uncertainty in the masses of the clus- 
ters. This result is essentially independent of the data we 
have used thus far and is not affe cted by varying A'cff. 

As can be seen in Figure 13 the cluster data pre- 
fer a value of tTg(Om/0.25)°''* that is lower than those 
preferred by SPT-hWMAP-hifo+BAO in the dn,/d\nk, 
free-1^, or free-A'eff models. Put another way, the con- 
straints on these parameters move closer to their stan- 
dard values when the cluster information is included. 
This effect is most significant for the constraint on the 
effective number of relativistic species, which moves to 
iVefi = 3.42 ±0.32 and is 1.2cr from the standard value of 
A'eff = 3.046. Similarly, the spectral running and primor- 
dial helium abundance are moved to within 1.4cr of their 
standard values when the cluster information is included. 
We conclude that models with a negative spectral run- 
ning, a high Yp, or a high N^s are disfavored by the clus- 
ter abundance data. The full parameter constraints for 
these models using SPT-HWMAP+i7o+BAO+Clusters 
are given in Table ro] 

We also note thatlower values of as are obtained if neu- 
trinos are allowed to have mass. For example, for a model 
in which A'efi and ^ rriv are allowed to be free, we find 
A^eff = 3.98 ±0.43, as = 0.803 ±0.056, and < 0.69 

eV (95% U.L.) using SPT±WMAP+i?o±BAO (com- 
pared to A^eff = 3.86 ± 0.42 and erg = 0.871 ± 0.033 if 
neutrinos are forced to be massless). 

5. CONCLUSION 

We have presented a new measurement of the damp- 
ing tail of the CMB power spectrum using data from the 
South Pole Telescope. This measurement builds upon 
earlier measurernents of the d amping tail by ACBAR 
(iReichardt et al. 120091), QU aD (iBrown et al.|2 009l|Fried-| 
manet al.||2009| ) and ACT (Pas et al. 2Ulla|- 'i'hg"5FT 
power spectruni uses 150 GHz data and spans the mul- 
tipole range 650 < I < 3000, where it is dominated by 
primary CMB anisotropy. We combine this spectrum 
with data from WMAP7 to constrain cosmological mod- 



els. We find that the SPT and WMAP7 spectra are con- 
sistent with each other, and that when combined they 
are well fit by a spatially flat, ACDM cosmology. 

The addition of the SPT data provides modest im- 
provements to the constraints on the standard six- 
parameter model relative to using WMAP alone. One 
notable improvement is that SPT+WMAP7 measure the 
scalar spectral index to be Ug = 0.9663 ± 0.0112, which 
disfavors the Harrison-Zel'dovich-Peebles index {ug = 1) 
at 3. Oct using only CMB data. When lo w-redshift mea- 



surements of the Hubble constant (Ricss et al. 
and the BAO feature fPercival et al.|201Ul are inc 



20111 
uded, 

the constraint on the scalar spectral index improves to 
Us = 0.9668 ± 0.0093, a 3.6ct rejection of Ug — 1. 

We consider a number of extensions beyond this base- 
line model. First we consider a model in which the am- 
plitude of gravitational lensing on the CMB is allowed 
to vary freely, and find that the SPT+WMAP data de- 
tect, at ~5ct, the effect of gravitational lensing, and that 
the amplitude is consistent with the ACDM cosmologi- 
cal model. Parametrized in terms of a rescaling of the 
lensing potential power spectrum {Cf'^ — >■ Ai^Cf'^), the 
lensing amplitude is A^,*'^ = 0.94 ± 0.15. 

We consider a model in which the power from ten- 
sor fluctuations is allowed to vary freely. We constrain 
the tensor-to-scalar ratio to be r < 0.21 (95% CL) 
using SPT+WMAP7, and r < 0.17 (95% CL) using 
SPT+WMAP7+i/o+BAO. 

We consider a model in which the scalar spectral index 
Ug is allowed to vary or "run" as function of wavenumber. 
We constrain the spectral running to be dug/dhik = 
-0.024 ± 0.013 using SPT+WMAP7. 

We consider a model in which the primordial helium 
abundance, typically a function of standard BBN theory, 
is allowed to vary freely. That is, we measure the effect 
of helium due solely to its effect on the CMB damping 
tail. We strongly detect the effect of helium on the CMB; 
a model with no helium is rejected at 7.7ct. When the 
primordial helium abundance is allowed to vary freely, 
we find Yp = 0.296 ± 0.030 using SPT+WMAP7. 

Finally, we consider a model in which the effec- 
tive number of relativistic species in the early uni- 
verse is allowed to vary freely. Normally this is the 
number of neutrinos, three, plus a small correction 
due to electron-positron energy injection, resulting in 
^sbindard ^ 3 Q4g^ ^siug SPT+WMAP7 wc strongly 
detect the effect of neutrinos on the CMB; a model with 
no neutrinos is rejected at 7.5ct. When A^eff is allowed 
to vary freely, we find A^off — 3.85 ± 0.62, while using 
SPT+WMAP7+i/o+BAO we find Aeff = 3.86 ± 0.42. 

Three of these model extensions — spectral running, 
free helium, and free A'eff — show a mild, ~1.7ct prefer- 
ence for non-standard models. We find that such models 
are disfavored by the value of CTg inferre d from the abun- 
dance of low-redshift galaxy clusters (Vikhlinin et al. 
2009 1 . The constraints on these parameters move closer 
to their standard values when the cluster information 
is included. Using SPT±WMAP7+i?o±BAO±Clusters, 
the constraints are dus/dlnk = — 0.017 ± 0.012, Yp = 
0.288 ± 0.029, and A^eff = 3.42 ± 0.32. 

The SPT data presented here cover 790 square de- 
grees. The full SPT-SZ survey, which is expected to be 
completed by the end of 2011, will cover approximately 
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Fig. 12. — The two-dimensional marginalized constraint on the primordial helium abundance Yp and the effective number of relativistic 
species A^^ff for a model in which both parameters are free. The two-dimensional contours show the 68% and 95% confidence intervals. 
The relation between the two quantities in standard BBN theo ry i s shown by the dashed line, with the point {N^g = 3.046, Yp = 0.2478) 
shown by the square. The constraint on N^g shown in Figure |11| is essentially a cut through this likelihood along the BBN curve, while 
the constraint on Yp shown in Figure [To| is a cut along N^ff = 3.1)46. 
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Fig. 13. — The two-dimensional marginalized constraints on spectral running, primordial helium, or the effectiv e number of relativist 
species versus the combination eg (nM/0.25)'''*^, which is well constrained by the cluster abundance measurement of |Vikhlinin et al. (2009 
Each panel corresponds to a distinct Markov chain. "CMB" corresponds to SPT-fWMAP7. The two-dimensional contours snow tlie bb' 
and 95% confidence intervals. The constraint on (Tg(f2jv//0.25)'' '* from the clusters and the corresponding la uncertainties are shown 
by the vertical lines. The standard values of the spectral running, primordial helium, and the effective number of relativistic species are 
shown by the dotted horizontal lines. Adding the cluster abundance information moves the constraints on these parameters closer to their 
standard values. 
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2500 square degrees. With 150 GHz data of the quahty 
used here and with additional data at 90 and 220 GHz, 
a power spectrum analysis of the full SPT survey should 
be at least 1.7 times more sensitive than that presented 
here. 
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